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ABSTRACT 

Context. The nearby high-mass star binary system ^^OriC is the brightest and most massive of the Trapezium OB stars at the core of the 
Orion Nebula Cluster, and it represents a perfect laboratory to determine the fundamental parameters of young hot stars and to constrain the 
distance of the Orion Trapezium Cluster 

Aims. By tracing the orbital motion of the ^^Ori C components, we aim to refine the dynamical orbit of this important binary system. 
Methods. Between January 2007 and March 2008, we observed 6/^OriC with VLTI/AMBER near-infrared {H- and /^T-band) long-basehne 
interferometry, as well as with bispectrum speckle interferometry with the ESO 3.6 m and the BTA 6 m telescopes {B'- and V'-band). 
Combining AMBER data taken with three diff'erent 3-telescope array configurations, we reconstructed the first VLTI/AMBER closure-phase 
aperture synthesis image, showing the ^^OriC system with a resolution of ~ 2 mas. To extract the astrometric data from our spectrally 
dispersed AMBER data, we employed a new algorithm, which fits the wavelength-diff'erential visibility and closure phase modulations along 
the H- and A'-band and is insensitive to calibration errors induced, for instance, by changing atmospheric conditions. 

Results. Our new astrometric measurements show that the companion has nearly completed one orbital revolution since its discovery in 1997. 
The derived orbital elements imply a short-period (P ^ 11.3 yr) and high-eccentricity orbit (e ^ 0.6) with periastron passage around 2002.6. 
The new orbit is consistent with recently published radial velocity measurements, from which we can also derive the first direct constraints on 
the mass ratio of the binary components. We employ various methods to derive the system mass (Msystem = 44 + 7 M©) and the dynamical 
distance (d = 410 + 20 pc), which is in remarkably good agreement with recently published trigonometric parallax measurements obtained 
with radio interferometry. 
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1. Introduction 

The Orion Nebula (M42) is one of the closest (d ~ 400 - 450 
pc) and most prominent star-forming regions (see O'Dell 2001 
for a review). It contains a massive cluster of very young 
(~ 1 X 10^ yr) stars (cf. Herbig & Terndrup 1986; McCaughrean 
& StauflTer 1994; Hillenbrand 1997), which is known as the 
Orion Nebula Cluster (ONC). Due to its relatively close dis- 
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* Based on observations made with ESO telescopes at the La 
Silla Paranal Observatory under the OT and VISA-MPG GTO pro- 
gramme IDs 078.C-0360(A), 080.C-0541(A,B,C,D), 080.D-0225(B), 
and 080.C-0388(A). 



tance and its favorable celestial location (which makes it ac- 
cessible to observatories on the northern and southern hemi- 
spheres), the ONC is probably the best investigated young clus- 
ter in the whole sky and has been observed at virtually every 
wavelength. It is a perfect laboratory for observations of young 
stellar objects over the full mass range, from very low-mass 
brown dwarfs to massive 0-type stars. 

The brightest star in the cluster is the massive 07-05.5 
type ^ star ^^OriC, which is known to be a close visual bi- 

^ Simon-Diaz et al. (2006) determined the stellar eff'ective temper- 
ature of e^OriC to be T^q = 39000 + 1 000 K with logg = 4.1 dex. 
In recent stellar atmosphere models (Martins et al. 2005), this cor- 
responds to a spectral type of 05.5-06, while the same parameters 
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nary system. After the initial discovery of the companion at a 
separation p of Of! 033 (33 mas, corresponding to about 15 AU) 
with near-infrared (NIR) bispectrum speckle interferometry by 
Weigelt et al. (1999), Schertl et al. (2003) have presented fur- 
ther observations and reported the first detection of orbital mo- 
tion. Kraus et al. (2007) presented the first speckle observa- 
tions at visual wavelengths, the first NIR long-baseline inter- 
ferometric observations of 6^ Ori C using the IOTA interferom- 
eter, and produced an aperture-synthesis image of the system. 
They also performed a joint analysis of all existing interfero- 
metric measurements that covered a period of more than 9 years 
and clearly revealed orbital motion. After reaching a maximum 
value of 42 mas in 1999, the separation of the system steadily 
decreased to 13 mas in 2005. Detailed modeling of these data 
yielded a preliminary orbit solution with a high eccentricity 
(e ~ 0.91) and a period of 10.9 yrs. According to this solu- 
tion, the periastron passage should have occurred around July 
2007 with a closest separation of less than 2 AU. Patience et al. 
(2008) recently presented additional interferometric observa- 
tions of ^^OriC obtained with NPOI at visual wavelengths. 
Extending the orbital coverage by about 1 .2 yrs, they measured 
a companion position which deviates ~ 4 mas from the posi- 
tion predicted by the orbital solution of Kraus et al. (2007) and 
concluded that the orbit has a considerably lower eccentricity 
(e ~ 0.16) and a longer period (~ 26 yrs). This shows a clear 
need for further interferometric observations that will clarify 
this apparent discrepancy. A reliable and accurate orbit solution 
will also provide unique constraints to the basic stellar param- 
eters and, most importantly, will yield the masses of the stars 
and the distance. 

The orbital motion of the companion might also be respon- 
sible for the radial velocity variations reported by Stahl et al. 
(1996, 2008) and others. Due to the incomplete coverage and 
the significant scatter in the derived velocities, it is not yet pos- 
sible to derive the orbital elements of the spectroscopic orbit 
(Stahl et al. 2008). 

In spite of its importance and the huge number of studies, 
the distance to the ONC is, even after decades of investiga- 
tion, still not well known and an issue of ongoing discussion. 
The "canonical" value of 450 pc (Herbig & Terndrup 1986; 
Genzel & Stutzki 1989) that was widely used during the last 
two decades was recently challenged by some studies find- 
ing significantly smaller values. Numerous new distance de- 
terminations (e.g. Stassun et al. 2006; Jeff"ries 2007; Mayne 
& Naylor 2008), also including the first direct trigonometric 
parallax measurements on radio sources in the ONC by Hirota 
et al. (2007, 437 ± 19 pc), Sandstrom et al. (2007, 389!^| pc), 
and Menten et al. (2007, 414 + 7 pc), yielded distances mainly 
in the range between 390 pc and 440 pc. The ~ 13% diflfer- 
ence in these distance estimates, which to some extent might 
also include an intrinsic distance spread of the studied stars, 
leads to an ~ 30% uncertainty in the derived stellar luminosities 
and correspondingly aff'ects any age and mass estimates for the 
stars. This is a serious limitation to the usefulness of the excep- 
tionally well-studied young stellar population; e.g., for the cali- 



correspond to a later spectral type of 07.5-08 using earlier models 
(e.g. Vacca et al. 1996). 



Table 2. AMBER calibrator stars and their characteristics, including 
uniform disk (UD) diameters. 



Star 


V 


H 


K 


Spectral 


Adopted UD diameter 










Type 


[mas] 


HD 33833 


5.9 


3.9 


3.8 


G7III 


0.83 + 0.06" 


HD 37128 


1.7 


2.4 


2.3 


BOI 


0.86 + 0.16^ 


HD 41547 


5.9 


5.1 


5.0 


F4V 


0.41 + 0.03^ 


HD 43023 


5.8 


3.7 


3.5 


G8III 


0.98 + 0.07^ 


HD 50281 


6.6 


4.3 


4.1 


K3V 


0.77 + 0.10^ 



Notes - UD diameter computed with ASPRO 

(http : //vivivi . jrmtic . f r/aspro_page . htm) . 

^ UD diameter taken from Mozurkewich et al. (1991). 

^ UD diameter taken from getCal tool 
(http : //mscweb . ipac . caltech . edu/gcWeb/) . 



bration of pre-main-sequence evolutionary models or in the in- 
vestigation of the spatial relationship between the young stellar 
groups in the wider region of the Orion association. Therefore, 
trigonometric parallax measurements of stars in the core of 
the Trapezium OB star cluster (e.g. on GMR 12=^^0ri A2, for 
which Menten et al. 2007 measured 418 + 9 pc) are highly de- 
sirable. An alternative way to obtain such distance estimates 
is through astrometric measurements on close binary systems 
such as ^^Ori C, yielding the dynamical parallax of the system. 

2. Observations and data reduction 

2.1. Bispectrum speckle interferometry 

Bispectrum speckle interferometry is a powerful technique 
to overcome atmospheric perturbations and to reach the 
diffraction-limited resolution of ground-based telescopes at vi- 
sual wavelengths. After monitoring the orbital evolution of the 
system between 1997 and 2004 at wavelengths between 2. 1 fim 
and 545 nm (p = 43 . . . 24 mas; Weigelt et al. 1999; Schertl 
et al. 2003; Kraus et al. 2007) using the Russian Big Telescope 
Alt-azimuthal (BTA) 6.0 m telescope, ^^OriC became unre- 
solvable with 6 m-class telescopes in 2006. In November 2007 
(p = 20 mas) and January 2008, the system could again be re- 
solved with the BTA 6 m telescope using a V medium-band 
filters {Ac = 550 nm, A A = 20 nm). For these observations, a 
1280x1024 pixel CCD with a multi-alkali S25 intensifier pho- 
tocathode was used. 

In January 2008, additional speckle measurements using V 
(Ac = 550 nm, AA = 30 nm) and B' (Ac = 440 nm, AA=16 nm) 
medium-band filter were obtained with the ESO 3.6 m tele- 
scope. The speckle observations were carried out in the course 
of ESO open-time progranmie 080.C-0388(A) using our visi- 
tor speckle camera, which employs a Marconi/EEV electron- 
multiplying CCD. For the speckle observations, we recorded 
interferograms of ^^Ori C and of nearby unresolved stars in or- 
der to compensate for the atmospheric speckle transfer func- 
tion. The calibrator stars, the number of recorded interfero- 
grams and the detector integration times (DITs) are listed in 
Table 1. The modulus of the Fourier transform of the ob- 
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Table 1. Observation log for our new bispectrum speckle and long-baseline interferometric observations. 



Instrument 


Date 


UT 


Telescope 


Spectral 


DIT 


No. Interferograms 


Calibrator(s) 




[UT] 




Triplet 


Mode 




Target/Calibrator 




B 1 A m/5pecKle 


zUU/ JNov 1j 






jjUnm/zU nm 


20 ms 


194U/194U 


u Uri D 


b5(j J.om/^pecKle 


zUUo Jan lU 






440 nm/16nm 


10 ms 


loUUU/lUUUU 


JO (Jri 


J.om/^pecKle 


zUUo Jan lU 






jjUnm/JUnm 


10 ms 


lUUUU/oUUU 


Jo Uri 


B 1 A D m/5pecKle 


zUUo Jan zb 






jjUnm/zU nm 


20 ms 


1 jUU/zUUU 


u (Jri D 


VLTl/ AMBER 


2007 Jan 08 


04:36, 06:32 


T TT' 1 T TTO T TT A 

UT1-UT3-UT4 


LR-^ 


26 ms 


12000/10000 


T TT~V A t C A ^ 

HD 41547 


VLTl/AMBbR 


2007 Dec 03 


07:19 


AO-DO-HO 


LR-//^ 


26 ms 


5000/10000 


HD 33833 


VLTI/ AMBER 


2007 Dec 03 


07:34, 08:44 


AO-DO-HO 


LR-//^ 


50 ms 


7500/7500 


T TTA O O O O O 

HD 33833 


VLTl/ AMBER 


2007 Dec 03 


r\'~i A o r\o c A 

07:48, 08:54 


AO-DO-HO 


LR-//^ 


100 ms 


2000/750 


T TTA O O O O O 

HD 33833 


\ Tl TPT / A A TT> "Cn 

VLTl/ AMBER 


2007 Dec 05 


06:05 


AO-KO-Gl 




26 ms 


C AA A ICf\f\f\ 

5000/5000 


T TTA O O O O O 

HD 33833 


VLTl/ AMBER 


2007 Dec 05 


05:49, 07:28 


A f\ 'XT' f\ 1 

AO-KO-Gl 


LR-HK 


50 ms 


10000/15000 


T TTA O O O O O 

HD 33833 


VLTl/AMBER 


2007 Dec 05 


06:17,07:46 


AO-KO-Gl 


LR-HK 


100 ms 


3500/2500 


HD 33833 


VLTl/AMBER 


2008 Feb 22 


03:44 


AO-DO-HO 


LR-HK 


50 ms 


5000/10000 


HD 37128, HD 50281 


VLTl/AMBER 


2008 Feb 22 


04:21 


AO-DO-HO 


LR-HK 


100 ms 


2500/3500 


HD 37128, HD 50281 


VLTl/AMBER 


2008 Feb 24 


03:11 


AO-KO-Gl 


LR-HK 


50 ms 


10000/5000 


HD 37128 


VLTl/AMBER 


2008 Feb 24 


02:37 


AO-KO-Gl 


LR-HK 


100 ms 


6000/4000 


HD 37128 


VLTl/AMBER 


2008 Mar 03 


02:31 


DO-HO-Gl 


LR-HK 


50 ms 


5000/5000 


HD 43023 


VLTl/AMBER 


2008 Mar 03 


02:52 


DO-HO-Gl 


LR-HK 


100 ms 


1500/1500 


HD 43023 



ject (visibility) was obtained with the speckle interferometry 
method (Labeyrie 1970). For image reconstruction we used 
the bispectrum speckle interferometry method (Weigelt 1977, 
Weigelt & Wirnitzer 1983, Lohmann et al. 1983, Hofmann & 
Weigelt 1986). 

Besides providing an independent astrometric measure- 
ment, our bispectrum speckle interferometric measurement 
is of special importance as it allows us to solve the 180°- 
ambiguity, which is inherent to long-baseline interferometric 
investigations which do not include the closure phase (e.g. 
Patience et al. 2008) or for which the instrumental closure 
phase sign has not yet been calibrated (as for VLTl/AMBER). 
Solving this ambiguity is essential for deriving the orbit of 
the system. Therefore, we paid special attention while deriv- 
ing the field orientation of our speckle images, using calibra- 
tion measurements which were taken with the same instrument 
setup during the same night as the ^^OriC observations on 
^^Ori A, ^^OriB and the well-studied multi-component object 
T] Carinae. For the January 2008 observations, a position angle 
calibration with an accuracy of ~ 0.3° was done using ^'-band 
observations covering the Trapezium stars ^^Ori A, B, and E, 
followed by the V and observations of ^^OriB, which is 
a binary system with a separation of -T', where the fainter 
component is a close binary with a separtion of -O'.'IS. This 
allows us to unambiguously determine that in January 2008, 
the fainter component (C2) was located to the southwest of the 
primary star (Fig. 1 and A.l). Performing a detailed re-analysis 
on all Speckle data taken by our group between 1997 and 2005, 
we found a 180° -calibration problem which affected the po- 
sition angle of the Speckle measurements at epochs 2003.8, 
2003.9254, 2003.928, and 2004.8216, as published in Kraus 
et al. 2007. Since these Speckle measurements were used for 
the calibration of the IOTA and NPOI long-baseline interfero- 
metric observations by Kraus et al. (2007) and Patience et al. 
(2008), a revision of the preliminary orbital solutions presented 
in these papers is required (as presented in Sect. 5.1). In Tab. 3, 
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Fig. 2. wv-plane coverage obtained at four epochs on ^^OriC with 
VLTl/AMBER. The radial extension of the wv-tracks reflects the spec- 
tral coverage of our AMBER interferometric data, covering either the 
K-hmd (LR-K mode, 2007 Jan 08) or the H- and K-bmd (LR-HK 
mode, 2007 Dec 03-05, 2008 Feb 22-24, 2008 Mar 03). 



we list all available astrometric data, taking the quadrant cor- 
rection into account. 

2.2. VLTl/AMBER spectro-interferometry 

AMBER is the NIR interferometric beam-combiner instrument 
(Petrov et al. 2007) of the Very Large Telescope Interferometer 
(VLTI), which is located on Cerro Paranal/Chile and operated 
by ESO. Eor the beam combination, either three 8.2 m unit tele- 
scopes (UTs) or three of the allocatable 1.8 m auxiliary tele- 
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Fig. 1. Selection of interferometric images of the ^^OriC system obtained by our group between 1997 and 2008, revealing the orbital motion 
of the companion. The images were reconstructed either from H-, 7-, V-, or 5' -band speckle interferograms recorded at the BTA 6 m or the 
ESO 3.6 m telescopes, or from //-band Michelson interferograms recorded with the IOTA 3-telescope interferometer (see Kraus et al. 2007 
for details). For epoch 2007.9, we show the aperture synthesis image which we reconstructed from VLTI/AMBER data using the procedure 
described in Sect. 4. For each image, 10% intensity level contours are shown. In the speckle images, the fitted component positions are marked 
with an asterisk. 



scopes (ATs) can be used. One outstanding feature of AMBER 
is its spectral capability, allowing one to observe several spec- 
tral bands with low spectral resolution simultaneously (LR- 
JHK mode, R = A/AA = 35, covering the /-band around 
1.2 yum, the //-band around 1.6 yum, and the ^-band around 
2.2 yum). The wide wavelength coverage offered by this spectral 
mode allows us to fit the binary separation vector with a high 
accuracy and observing efficiency (as discussed in Sect. 3). 



The VLTI/AMBER data was recorded in the course of three 
ESO programmes, yielding the wv-plane coverage shown in 
Fig. 2. For each science observation, we also recorded inter- 
ferograms for at least one interferometric calibrator star (see 
Tab. 2 for the characteristics of our calibrator stars), which al- 
lows us to calibrate the visibilities and the closure phases for 
instrumental as well as atmospheric effects. 
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For AMBER data reduction, we employed the amdlibl- 
data reduction software ^ (release 2.1), which is based on the 
P2VM algorithm (Tatulli et al. 2007). For the LR-if ^ data, we 
applied the data selection criteria outlined in Kraus et al. (2008) 
in order to remove frames which were either degraded by atmo- 
spheric effects or were recorded significantly off'set from zero 
optical path delay. 

For the accuracy achievable with our binary model fits, 
the absolute calibration of the wavelength scale is of spe- 
cial importance. Therefore, we carefully refined the wave- 
length calibration using the telluric gaps between the ob- 
served spectral bands. Using this procedure (which is described 
in Appendix A), we reach a calibration accuracy of about 
0.03 fjm, which is still the dominating limiting factor on the 
total achievable astrometric accuracy (~ 2%). 

While the target/calibrator observations from January 2007 
were taken under good and stable atmospheric conditions, 
some of the observations in December 2007 and February 2008 
suff'er from strongly variable seeing conditions and short atmo- 
spheric coherence times. It is known that short coherence times 
can result in a decrease of the measured fringe contrast, which 
might not be completely calibrated out using calibrator mea- 
surements. This eff'ect can impose errors on the absolute cali- 
bration, which are expected to increase with longer DIT and to- 
wards shorter wavelengths. To illustrate this eff'ect, in Fig. 4 we 
plot the calibrated visibilities measured during our December 
2007 observation campaign and compare them to the cosine 
visibility modulation of a binary source. In particular, in the 
^-band the resulting calibration errors can be on the order of 
20%. 

However, since all spectral channels of an AMBER in- 
terferogram are recorded at the same time, the wavelength- 
differential observables (in particular, the differential visi- 
bility Sy{Xj) are practically insensitive to this degradation. 
Therefore, in Sect. 3 we employ a fitting algorithm in which 
we fit only differential visibilities and closure phases in order 
to determine the 6^ Ori C binary separation vector. 

A fundamental problem one encounters when interpreting 
VLTI/ AMBER 3 -telescope data concerns the lack of a calibra- 
tion measurement for the closure phase sign, resulting in a 180° 
uncertainty on the position angle of binary star observations. 
Since our ^^Ori C VLTI/ AMBER observations from December 
2007 and February 2008 bracket the bispectrum speckle mea- 
surement from January 2008, our data set allows us to un- 
ambigously define the closure phase sign for these observa- 
tions, providing a direct calibration of the 180° uncertainty for 
VLTI/AMBER for the first time. 

3. Modeling 

In order to derive the binary separation p and position angle ^ 
from the speckle data, we used the same algorithm as in our 
earlier studies on this system (Weigelt et al. 1999; Schertl et al. 

^ The amdlibl software package is available from the website 
http : //www . jmmc . fr/data_processing_amber . htm 

^ Following convention, we measure the position angle (PA) East of 
North. 
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10 20 30 40 50 60 70 80 



x[10^:V] 

Fig. 4. AMBER visibilities measured on 2007 Dec 03 and 
2007 Dec 05 plotted versus projected distance in the wv-plane (where 
the projection was perpendicular to the fitted binary separation vec- 
tor = 241.2°, i.e. x = wcos(0 + 90°) - vsin(0 + 90°)). The solid 
red line shows the theoretical cosine visibility profile for a binary star 
with separation 19.07 mas and intensity ratio 0.30. As indicated by 
the strong visibility off'sets, which particularly occur at high spatial 
frequencies, the absolute calibration is sometimes rather poor, reflect- 
ing the changing atmospheric conditions during these nights. It can 
also be seen that this calibration bias is particularly important for long 
DITs (50 ms, 100 ms), while it is nearly negligible for short DITs 
(26 ms, grey & black points). As expected, the spectral dependence of 
the visibility is not aff'ected by these calibration uncertainties. 

2003; Kraus et al. 2007), fitting a cosine function directly to 
the 2-D speckle visibilities. For illustration, in the Appendix 
(Fig. A.l) we show power spectra and Fourier spectra deter- 
mined from three independent data sets obtained during our 
observing run with the ESO 3.6 m telescope and compare them 
with the model power and Fourier spectra corresponding to a 
binary star. 

For the modeling of the AMBER data, we employed an 
optimized algorithm which fits closure phases (CP) and dif- 
ferential visibilities dV. Each AMBER LR-//^ measurement 
records 16 ^-band plus 11 //-band spectral channels in the 
wavelength range from 1.51 to 2.55 yum, covering about 40% 
of the object Fourier spectrum in radial direction (see wv-plane 
tracks in Fig. 2). Strongly resolved objects (such as binaries 
with a separation of p > AA/B, where AA is the recorded spec- 
tral window and B' is the projected baseline length) can already 
show significant visibility modulation over this range of spatial 
frequencies (see illustration in Fig. 3). For instance, our ^^Ori C 
AMBER measurements from December 2007 probe spatial fre- 
quencies out to the eighth lobe of the cosine binary visibility 
modulation (Fig. 4) and the wavelength-diff'erential visibilities 
recorded in a single AMBER //-//T-band measurement sample 
up to 3 visibility cycles. This wavelength-diff'erential visibil- 
ity modulation already provides all information required to de- 
termine the binary separation and orientation. In this way, the 
absolute calibration of the visibility, which is subject to many 
adverse atmospheric eff'ects, becomes dispensable. 

To determine the diff'erential visibilities SV as used in our 
fitting algorithm, we first follow the standard data reduction 
and calibration procedure in order to correct for wavelength- 
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Fig. 3. Illustration of the basic principle of fitting binary parameters using wavelength-differential interferometric observables. Left: The 
VLTI/AMBER 3-telescope interferometer measures the interferometric observables in various spectral channels from 1.4 to 2.4 fim (as in- 
dicated by the color of the dots) and towards three different position angles, probing different regions in the two-dimensional Fourier spectrum 
of the source brightness distribution. The figure shows the visibility spectrum for a binary with p = 19.07 mas and = 241.2° (as inferred 
for ^^OriC on 2007 Dec 03) and the wv-sampling obtained with one of our AMBER observations on the same date. Middle and Right: The 
visibilities and phases show a wavelength-differential modulation which is independent of the absolute calibration (see Sect. 3) and which can 
be fitted to analytical models. 



dependent instrumental effects (using a calibrator measurement 
taken during the same night). Then, we remove the absolute 
calibration by subtracting the average visibility of the consid- 
ered spectral window 

SV(A) = V(A) - <y(^))band, (1) 

yielding the differential visibility SV, which we compute sepa- 
rately for each spectral band (H and ^-band). 

A similar approach is applied to the model visibilities 
Vmodei(^) (which wc computc using equations 7, 11, and 12 
from Kraus et al. 2005): 

SVjnodcl(A) = a (Vmodel W) " b) , (2) 

where a and b are adjusted to minimize the residuals between 
the model visibilities and the measured wavelength-differential 
visibilities dV{X) before a Levenberg-Marquardt least-square 
fitting algorithm is used to determine the best-fit model param- 
eters. The measured closure phases are fit simultaneously with 
the differential visibilities, weighting each data point accord- 
ing to the determined error bars (see equations 8-10 in Kraus 
et al. 2005). Parameters in our binary star model are the binary 
separation p, the position angle 0, the intensity ratio FqiIFqi, 
and the uniform disk diameter of the components. Since our 
earlier speckle measurements indicated that i^c2/^ci is practi- 
cally constant in the NIR wavelength range (see Kraus et al. 
2007), we assume that the intensity ratio does not vary over 
the H- and ^-bands. For the model fitting, we fix the apparent 



diameters of the two stars to 0.22 mas for the primary star (cor- 
responding to 10.6 Rq) at a distance of 450 pc) and 0. 15 mas for 
the companion (7.2 7?q, using the spectral type determined by 
Kraus et al. 2007 and the stellar evolutionary models of Martins 
et al. 2005). Our AMBER LR measurements and best- fit mod- 
els are shown in Fig. 5 and 6. The derived astrometric data are 
given in Tab. 3. 



4. VLTI/AMBER aperture synthesis imaging 

Currently, the most commonly applied procedure to extract sci- 
entific information from optical interferometric data is model 
fitting, which requires the selection of an appropriate (geomet- 
rical or astrophysically motivated) model, whose parameters 
are then adjusted to fit the interferometric observables. Since 
this model fitting approach requires a priori knowledge about 
the source structure, it might, in some cases, not be applicable 
or might lead to biased results. Therefore, it is highly desirable 
that optical interferometers such as VLTI/AMBER have the ca- 
pability of recovering the source brightness distribution free of 
any assumptions. The aim of this section is to apply state-of- 
the-art aperture synthesis imaging techniques for the first time 
to real VLTI/AMBER data in order to independently confirm 
the scientific results obtained for ^^OriC in the last section, 
and, simultaneously, to demonstrate the imaging capabilities of 
VLTI/AMBER on a relatively well-studied astrophysical target 
with limited intrinsic complexity. 
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Fig. 5. Differential visibilities and closure phases measured with VLTI/AMBER in January 2007 (left), February 2008 (middle), and March 2008 
(right) on ^^Ori C. The solid lines show the best-fit model corresponding to the binary parameters given in Tab. 3. 



In order to obtain the wv-coverage required for aperture syn- 
thesis imaging, we combined the ^^OriC AMBER data sets 
taken between December 2007 and March 2008 on three differ- 
ent 3 -telescope array configurations and at several hour angles. 
Assuming that the source morphology does not change signif- 
icantly over the ^-band (which seems well justified based on 
our earlier measurements of the wavelength-dependent binary 
flux ratio; Kraus et al. 2007), we make use of AMBER 's spec- 
tral coverage, yielding radial tracks in the wv-plane (Fig. 7, left). 
Since the /sT-band visibilities provide a more reliable absolute 
calibration (see Sect. 2.2), we did not use the //-band data for 
image reconstruction and rejected also one measurement taken 
under particularly poor and variable conditions (2007 Dec 05, 
UT 07:46). In order to correct for the binary orbital motion 
over the 3-month period (Ap = 2.2 mas, A0 = 6.6°), we 
apply a rotation-compensating coordinate transformation of 
the wv-plane (see Kraus et al. 2005 for a description of this 
procedure) using the astrometric data given in Tab. 3. Then, 
we employed our Building Block Mapping software, which 
is based on the algorithm described by Hofmann & Weigelt 
(1993). Starting from an initial single ^-function, this algo- 



rithm adds components to a model image in order to mini- 
mize the deviations between the measured bispectrum and the 
bispectrum of the model image. Finally, the image is con- 
volved with a clean beam of ~ 1.5 x 3 mas, reflecting the 
elongation of the sampled wv-plane. The resulting ^-band aper- 
ture synthesis image (Fig. 1 and Fig. 7, right) yields a direct, 
model-independent representation of our VLTI/AMBER inter- 
ferometric data. Measuring the companion position in the re- 
constructed image (p = 19.3 mas, = 241°, epoch 2007.9), 
we find good agreement with the results obtained with our 
model fitting approach using wavelength-diflTerential observ- 
ables (Sect. 3). The noise features within the image show an 
intensity amplitude below 2% of the peak brightness in the im- 
age. 

5. Results 

5. ) . Dynamical orbit of the 6^ Oh C binary system 

Since our new astrometric data extend the orbital coverage for 
the ^^OriC system by about 12 months and solves the calibra- 
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Fig. 6. Similar to Fig. 5, showing the AMBER data and best-fit model for our observations from December 2007 



tion problem described in Sect. 2.1, we can significantly im- 
prove the orbital solution of the system. To derive a refined 
orbital solution, we applied two independent orbit fitting ap- 
proaches. 

First, we used the method presented by Docobo (1985), 
which generates a family of Keplerian orbits, whose apparent 
orbits pass through three base points. These base points might 
be selected from the actual astrometric measurements or rep- 
resent observationally favored areas in the (p, 0, t) parameter 
space. From the generated class of possible solutions, the or- 
bit which best agrees with the measured separations and PAs is 
selected. We use the error bars of the individual measurements 
as weight. The orbital elements for the determined best-fit orbit 
solution are given in Tab. 4. 

In addition, we implemented a grid search algorithm which 
scans the parameter space of the dynamical elements eccen- 
tricity e, period P, and time of periastron passage Tq. As de- 
scribed by Hilditch (2001), at each grid point, the geometri- 
cal orbital elements a, /, fl, and oj can be determined by fit- 
ting the Thiele-Innes constants to the observational data. We 
scanned the parameter space between P = 9 . . . 16 yrs (in in- 



crements of 0.001 yrs), e = 0.0 ... 1.0 (in increments of 0.005), 
and To = 2000.0 . . . 2010.0 (in increments of 0.001 yrs) and 
determined the least- square distance between the N measured 
positions (p/, 0/) and the corresponding orbit positions (pj, ©J) 
for each orbit solution: 




Using the Xr value determined at each grid point, we built a 
Xr(P^ ^^Tq) data cube, which we projected to determine the lo- 
cal minimum and the associated uncertainty for each parame- 
ter. Fig. 9 shows the determined minima curves, from which 
we determine the best-fit orbital elements given in Tab. 4. Given 
that the orbital coverage has been substantially improved since 
our 2007 study, we do not have to impose a priori constraints 
on the distance/system mass in order to eliminate unphysical 
orbit solutions. In Fig. 8c we compare the orbital solutions de- 
termined with the Docobo (1985) algorithm and with the grid 
search algorithm to the available astrometric data. 
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Telescope 


Date 




Filter 


Flux ratio 


O 


P 


Ref. 












n 


[mas] 




BTA 6 m/Speckle 


1997, 


.784 


H 


0.26 + 0.02 


226.0 + 3 


33 + 2 


a 


BTA 6 m/Speckle 


1998, 


.838 


K' 


0.32 + 0.03 


222.0 + 5 


37 + 4 


a 


BTA 6 m/Speckle 


1999, 


.737 


J 


0.31 + 0.02 


214.0 + 2 


43 + 1 


b 


BTA 6 m/Speckle 


1999, 


.8189 


G' 


0.35 + 0.04 


213.5 + 2 


42 + 1 


c 


BTA 6 m/Speckle 


2000, 


.8734 


V 


0.35 + 0.03 


210.0 + 2 


40+1 


c 


BTA 6 m/Speckle 


2001, 


.184 


J 


0.29 + 0.02 


208.0 + 2 


38 + 1 


b 


BTA 6 m/Speckle 


2003, 


.8 


J 


0.30 + 0.02 


19.3 + 2 


29 + 2 


c 


BTA 6 m/Speckle 


2003, 


.9254 


v 




19.0 + 2 


29 + 2 


c 


BTA 6 m/Speckle 


2003, 


.928 






19.1 + 2 


29 + 2 


c 


BTA 6 m/Speckle 


2004, 


.8216 


V 


0.34 + 0.04 


10.5 + 4 


24 + 4 


c 


IOTA 


2005, 


.92055 


H 


0.28 + 0.03 


^ A ^ ^ A 

342.74 + 2 


13.55 + 0.5 


c 


NPOI 


2006, 


.1486 


V 




332.3 + 3.5 


11.80 + 1.11 


d 


VLTI/ AMBER 


2007, 


.0192 


K 


0.31 ± 0.03 


274.9 + 1 


11.04 + 0.5 




NPOI 


2007, 


.1425 


V 




268.1 + 5.2 


11 f\ A O 1 

11.94 + 0.31 


d 


NPOI 


2007, 


.1507 


V 




272.9 + 8.8 


12.13 + 1.58 


d 


NPOI 


2007, 


.1753 


V 




266.6 + 2.1 


12.17 + 0.37 


d 


NPOI 


2007, 


.2055 


V 




265.6 + 1.9 


12.28 + 0.41 


d 


NPOI 


2007, 


.2137 


V 




263.0 + 2.3 


12.14 + 0.43 


d 


BTA 6 m/Speckle 


2007, 


.9014 


Y 


0.29 + 0.02 


238.0 + 2 


19.8 + 2 




VLTI/AMBER 


2007, 


.9233 


H+K 


0.24 + 0.07 


241.2+1 


19.07 + 0.5 




ESO 3.6 m/Speckle 


2008, 


.0274 


B' 


0.29 + 0.02 


237.0 + 3 


19.7 + 3 




ESO 3.6 m/Speckle 


2008, 


.0274 


V 


0.31+0.02 


236.5 + 3 


19.6 + 3 




BTA 6 m/Speckle 


2008, 


.0712 


Y 


0.31 + 0.02 


236.2 + 2 


20.1 + 2 




VLTI/AMBER 


2008, 


.1479 


H+K 


0.23 + 0.09 


234.6 + 1 


21.17 + 0.5 




VLTI/AMBER 


2008, 


.1726 


H+K 


0.26 + 0.06 


236.4 + 1 


21.27 + 0.5 





Notes - The position angles given in col. 5 are measured East of North and were corrected for the 180° -calibration problem described in 
Sect. 2.1. 

References - (a) Weigelt et al. 1999, (b) Schertl et al. 2003, (c) Kraus et al. 2007, (d) Patience et al. 2008. 



5.2. Constraining tlie binary mass ratio 



Besides the 15.424-day period which is associated with the 
wind from the primary, long-term radial velocity variations 
were also found (Vitrichenko 2002; Stahl et al. 2008). Using 
a large data set covering more than 15 yrs of spectroscopic ob- 
servations (plus three archival measurements, which extend the 
coverage to more than 64 yrs), Stahl et al. (2008) showed that 
these variations are consistent with the orbital motion of a high- 
eccentricity binary system. Although the strong scatter within 
the radial velocity measurements prevents us from solving for 
the precise spectroscopic orbit, the combination of these data 
with our new orbital solution can be used to provide a first 
direct constraint on the mass ratio of the components in the 
^^OriC system. 

In Fig. 10 we plot the available radial velocity data as a 
function of orbital phase using the orbital period P and time 
for periastron passage Tq determined independently from our 
interferometric measurements (Sect. 5.1). Using the method 
from Pourbaix (1998), we compute the radial velocity varia- 
tions corresponding to our full set of orbital elements (Tab. 4) 
and perform a least-square fit between the measured (v/) and the 
predicted (v[) radial velocities in order to determine the mass- 



ratio between the binary components. As least-square measure 
for N measurements, we use 



1 



Vn 



(4) 



where Vq is the velocity of the center of mass of the system, 
which is given by Vq •= (v?) - (v^. Since the radial velocities 
were extracted from various spectral lines (C iv, He ii, and O iii) 
and it is known that these lines can show, with respect to each 
other, systematic velocity off'sets on the order of 2 - 3 km s~^ 
(Stahl et al. 2008), we used 3 km s~^ as minimum velocity error 
(Ty. in order to avoid overweighting individual measurements. 
By varying the mass-ratio between the components, we find 



^ Mc2 \ -^-^HUpc 



-1 



(5) 



i.e. ^(450 pc) = 0.21 +0.04 or ^(414 pc) = 0.23+0.05, which is 
slightly lower than the value we derived earlier by modeling the 
wavelength-dependent binary flux ratio of the ^^OriC system 
(q = 0.45 + 0.15, Kraus et al. 2007). 

For the radial velocity of the center of mass, we deter- 
mine 23.6 kms"\ which is in good agreement with the helio- 
centric velocity of the Orion Molecular Cloud (~ 28 kms"\ 
O'Dell 2001). This might indicate that the relative motion of 
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Fig. 7. Combining AMBER data obtained on three telescope configurations (Left: wv-coverage), we reconstructed an aperture synthesis image 
of the ^^Ori C system with an eff'ective resolution of ~ 2 mas (Right). For a detailed description, we refer to Sect. 4. 



a) Orbital period P 
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Fig. 9. Minimum curves as function of the dynamical orbital parameters P, e, Tq. These curves were used to determine the best-fit orbit 
solution and to evaluate the uncertainties on the individual parameters (see Tab. 4). 



the ^^Ori C system with respect to the parental cloud is smaller 
than previously assumed (O'Dell 2001; Stahl et al. 2008). 

5.3. Dynamical masses and parallaxes 

Kepler's third law (Msystem • = a^/P^) relates a and P with 
the product of the system mass Msystem •= ^ci + and the 
cube of the parallax tt. Therefore, our astrometric measurement 
of a^/P^ directly constrains certain areas in the (tt, Msystem) pa- 
rameter space, as shown by the shaded area in Fig. 13. Since the 
elements of the spectroscopic orbit are currently only weakly 
constrained, it is not possible to directly separate the system 
mass and the dynamical parallax. However, several indepen- 
dent methods exist to disentangle these important parameters 
using additional information, such as the flux ratio of the com- 
ponents, their mass ratio, and/or their stellar parameters. 

The stellar parameters of ^^OriCl have already been stud- 
ied extensively with spectroscopy, placing the efl'ective tem- 
perature between 37 000 and 40000 K (e.g. Rubin et al. 1991; 
Baldwin et al. 1991; Simon-Diaz et al. 2006; Polyakov & 
Tsivilev 2007). In particular, Simon-Diaz et al. (2006) included 
non-LTE and line-blanketing effects (reff,ci = 39 000+ 1000 K) 



for their detailed spectroscopic modeling, which makes their 
results compatible with the 0-star stellar parameter calibration 
by Martins et al. (2005). 

Since each of the above-mentioned stellar or observational 
parameters is associated with certain assumptions and uncer- 
tainties, it seems advisable to take a number of alternative ap- 
proaches for deriving the underlying physical parameters in 
order to yield some insight into the associated uncertainties. 
Therefore, we follow three alternative approaches: 

a) Based on the efl'ective temperature determination of the pri- 
mary and the binary flux ratio, one can give a reasonable 
mass range for each component and, thus, the system mass 
(green- shaded area in Fig. 13). 

b) One can estimate the mass of the primary from the stellar 
temperature and then derive the companion mass using the 
mass ratio constraints obtained from the radial velocities 
(Sect. 5.2), yielding the orange curve in Fig. 13. 

c) Baize & Romani (1946) presented a method which solves 
for the system mass and the dynamical parallax using a 
mass-luminosity relation (MLR), the bolometric correc- 
tions of the components and their extinction-corrected ap- 
parent magnitudes. To evaluate the influence of the MLRs 
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Table 4. Orbital elements as determined with the algorithm from 
Docobo (1985) and our grid search algorithm (described in Sect. 5.1). 
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Fig. 8. Comparison of our new orbit solutions with the available astro- 
metric data. Each position measurement is connected to the orbit pre- 
diction with an 0-C vector (arrows). Furthermore, we show the lines 
of nodes in grey. The plots are centered on the primary star. North is 
up and east is to the left. 
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Fig. 10. The radial velocities measured towards ^^OriC plotted as 
function of orbital phase (assuming the orbital elements determined 
with our grid search algorithm; see Tab. 4). Besides the extensive data 
set by Stahl et al. (2008), we included the radial velocities measured 
by Struve & Titus (1944), Conti (1972), and Morrell & Levato (1991). 
The red line shows the radial velocities corresponding to our orbital 
solution with a system velocity Vq - 23.6 kms~^ (blue dashed line) 
and the mass ratio ^(414 pc) = 0.23 as determined from our least- 
square fit (Sect. 5.2). 
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algorithm 


algorithm 
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[yrs] 
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CL 
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Notes - Besides the orbital elements, we give the mass ratio 
(Sect. 5.2), dynamical distance, and system mass (Sect. 5.3), derived 
from both set of orbit elements. The dynamical distance and system 
mass was determined using the method from Baize & Romani (1946, 
method c in Sect. 5.3) and three different MLRs. When assuming an- 
other distance d\ the dynamical system mass Mci + Mc2 must be 
scaled by a factor (d' /ddyn)^- The mass ratio MciIMqi was also com- 
puted for the distance Jdyn. but can be converted to any other distance 
using equation 5. 
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Fig. 11. curve for our fit of the binary mass ratio to the available 
radial velocity data (assuming a distance of 414 pc). We assume the 
orbital elements (Tab. 4) derived from our interferometric measure- 
ments. 

on the result, we used three different MLRs; namely, from 
Baize & Romani (1946), Heintz (1978), and Demircan & 
Kahraman (1991) (corresponding to the three blue data 
points in Fig. 13). 

For methods a) and b), we employ the stellar calibration 
by Martins et al. (2005), while for method cj, three earlier 
calibrations are used. In method a) and b), the stellar tem- 
perature of the primary component is fixed to the value by 
Simon-Diaz et al. (2006, reff,ci = 39 000 ± 1000 K), while 
in c) we scan a slightly wider range of temperature values 
(7^eff,ci = 37 000 . . . 40 000 K, T^^^ci = 30 000 . . . 33 000 K, 
using the bolometric correction by Bessell et al. 1998). Method 
a) requires the V-band flux ratio (Fci/Fci = 0.31 ± 0.02) 
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Fig. 12. Minimum curve as a function of a^/P^ = (^system) • tt^, 
through which the dynamical orbital parameters can be related to the 
mass sum and the dynamical parallax. For illustration, we give the 
dynamical distance in the upper axis assuming a total system mass of 
45 Mq. 
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Fig. 13. The derived a^/P^ value puts some direct constraints on the 
system mass Msystem •= Mqi + Mc2 and the dynamical distance. As 
described in Sect. 5.3, we employ three alternative methods to obtain 
constraints, which help to disentangle the system mass and the dis- 
tance (method a: dashed green-shaded area, method b: orange curve, 
method c: blue data points). The black-shaded area marks the a^/P^- 
constraints determined with our orbit grid search algorithm. With the 
red ellipse, we visually mark the area of best agreement between all 
constraints. 

and c) the extinction-corrected magnitude of the total system 
(V = 5.12 + 0.1, Ay = 1.74 ± 0.1, Hillenbrand 1997). 

Evidently, each method is associated with considerable un- 
certainties, which makes it very desirable for future obser- 
vations to improve not only the astrometric orbit, but to de- 
rive the accurate spectroscopic orbit of the system as well. 
Nevertheless, within their large uncertainties, the methods em- 
ployed cover a common area in parameter space, corresponding 



to a system mass of Msystem = 44 + 7 Mq and a dynamical dis- 
tance ofd = 410 ± 20 pc, as marked with the red ellipse in 
Fig. 13. Methods a) and b) yield systematically lower system 
masses than method c), reflecting the recent correction in the 
mass calibration scale in stellar evolutionary models (Martins 
et al. 2005). Assuming a distance of 414 pc, as determined by 
Menten et al. (2007), would yield a system mass of 46 Mq. 

5.4. Possible implications on the dynamical history of 
the 0^ Oh C system 

With an eccentricity of ~ 0.6, the orbit of ^^OriC is located 
on the upper end of the eccentricity distribution of low- as well 
as high-mass binary stars (Mathieu 1994; Mason et al. 1998), 
perhaps providing important information about the dynami- 
cal history of the system. Tan (2004, 2008) proposed that the 
Becklin-Neugebauer (BN) object, which is located 45'' north- 
west of the Trapezium stars, might be a runaway B star ejected 
from the ^^Ori C multiple system approximately 4 000 yrs ago. 
This scenario is based on proper motion measurements, which 
show that BN and ^^Ori C recoil roughly in opposite directions. 
Three-body interaction is a crucial part of this interpretation, 
and the high eccentricity of the 6^ Ori C orbit which we derive 
for this system might be a direct consequence of this dynami- 
cal interaction event. However, another study (Rodriguez et al. 
2005) also aimed to identify the multiple system from which 
BN was ejected and identified Source I as the likely progenitor 
system. Later, Gomez et al. (2005, 2008) added further evi- 
dence to this interpretation by identifying Source ^ as a poten- 
tial third member of the decayed system. Therefore, it is still 
unclear whether this scenario can explain the measured proper- 
ties of the Ori C orbit. 

As an alternative explanation for the high eccentricity of 
the 0^ Ori C system, Zinnecker & Yorke (2007) pointed out that 
such systems are predicted by star formation scenarios which 
include sub-Keplerian rotating disks or filament fragmentation 
(Kratter & Matzner 2006). 



6. Conclusions 

We have presented new bispectrum speckle (y'-/5'-band) and 
VLTI/AMBER (H-fK-band) interferometric observations of 
the Orion Trapezium star 6^ Ori C covering several epochs over 
a time period of about 14 months. 

From our long-baseline interferometric data, we have re- 
constructed the first model-independent VLTI/AMBER aper- 
ture synthesis image, depicting the 0^ Ori C system at a resolu- 
tion of ~ 2 mas and demonstrating the imaging capabilities of 
this unique facility. In order to extract accurate astrometric in- 
formation for all epochs, we have followed a new modeling ap- 
proach which is based on wavelength-diff'erential observables 
and which demonstrates the benefits of spectro-interferometry 
in terms of observing efficiency and robustness to poor observ- 
ing conditions. Furthermore, our ESO 3.6m and the BTA 6 m 
speckle observations allow us to solve 180° -ambiguity and to 
calibrate the closure phase sign of our VLTI/AMBER observa- 



Kraus et al: Tracing the high-eccentricity binary ^^Ori C through periastron passage 



13 



- 












Atmospheric Transmission (R=4000) ; 
Atmosplieric Transmission (R=35) -f 

/ \ AMBER P2VM spectrum-'' 
/ AMBER (RAW-SKY)/P2VM spectrum 






















































- 








.Ill 




ill 




























ji 





















































1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

;i[^m] 



Fig.A.2. To obtain a spectrum which is corrected for instrumental 
effects for the spectral calibration, we subtract from the RAW spec- 
trum the SKY spectrum and then divide by the P2VM spectrum 
(blue dashed curve). The resulting instrument-corrected spectrum (red 
curve) is compared with an atmospheric transmission spectrum (green 
curve: R = 4000; black curve: R = 35). 

tions, providing a potential reference for other AMBER obser- 
vations using closure phase information 

Our new 6^ Ori C astrometric data shows that since its dis- 
covery in 1997, the ^^OriC companion has nearly completed 
one orbital revolution and that the system has a high eccentric- 
ity (e ~ 0.6). Solving for the orbital elements, we determine a 
period of ~1 1.3 yrs, a semi-major axis of 44 mas, and a perias- 
tron passage around 2002.6. According to our orbital solutions, 
the physical separation between the components decreases to 
~ 7 AU (2.8 mas) during periastron passage, which might be 
too large to result in detectable signatures of wind- wind inter- 
action between the stellar winds. 

Using additional information about the stellar parameters 
and various stellar models, we estimate the total system mass to 
be 44 ± 7 Mq and the dynamical distance to be J = 410 ± 20 pc. 
In the coming years, the uncertainties on these parameters 
could be considerably reduced with new spectroscopic and as- 
trometric observations of this important binary system. 
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Appendix A: Spectral calibration of AMBER-LR 
data 

Employing wavelength-differential visibilities and phases for 
model fitting requires a precise knowledge of the central wave- 
length of the recorded spectral channels. Therefore, we per- 
formed a re-calibration of the wavelength axes of our data sets 
using atmospheric absorption features and found deviations of 

4 Since our data set might be useful for VLTI/AMBER 
users as reference data for the calibration of the closure 
phase sign, we provide our data on the following website: 
http : //www . mpifr . de/staf f /skraus/f iles/amber . htm 



up to ~ 0.1 yum compared to the standard wavelength table ap- 
plied to the data on Paranal. 

In the recorded /-///-//sT-band object raw spectra, the tel- 
luric features are not pronounced enough to be used for the 
calibration, which is mainly due to absorption by internal op- 
tical components. Therefore, besides the raw spectrum of the 
object exposures, we also extracted the spectrum from the sky 
and the P2VM calibration exposures. The P2VM calibration 
files (F2YM=Pixel-to-Visibility-Matrix; see TatulH et al. 2007) 
are recorded at the beginning of each observation block using a 
lamp located in AMBER's Calibration and Alignment Unit and 
provide a measure of the instrumental transmission. By sub- 
tracting the sky spectrum from the object spectrum and then 
dividing by the P2VM spectrum, we yield a spectrum which 
is corrected for most instrumental eff'ects and shows the tel- 
luric absorption features much more clearly (see Eig. A.2, red 
curve). These corrected spectra were compared with standard 
atmospheric transmission spectra provided by the Gemini ob- 
servatory^ (see Eig. A.2, green curve). We manually align the 
spectra using the pronounced gaps between the J-/H- and H- 
/^-band and find best agreement assuming a linear dispersion 
law A(i) oc 0.0328 • /, where / is the number of the spectral 
channel on the detector. 

The remaining wavelength calibration uncertainty of about 
1 spectral channel (corresponding to about 0.03 jim) is the 
dominant limiting factor on the achievable astrometric accu- 
racy (~ 2%). Therefore, for future observations, it seems highly 
desirable to implement an AMBER on-site spectral calibration 
device (e.g. NIR lasers with well-known frequencies), fascil- 
itating an absolute spectral calibration in the standard instru- 
ment calibration procedure, pushing AMBER to its full astro- 
metric accuracy. 
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